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We consider the critical spin-spin correlation function of the 2D Ising model with a line defect 
which strength is an arbitrary function of position. By using path-integral techniques in the con- 
tinuum description of this model in terms of fermion fields, we obtain an analytical expression for 
s ! ' the correlator as functional of the position dependent coupling. Thus, our result provides one of the 

| few analytical examples that allows to illustrate the transit of a magnetic system from scaling to 

q_) . non-scaling behavior in a critical regime. We also show that the non-scaling behavior obtained for 

^"*^ ■ the spin correlator along a non-uniformly altered line of an Ising model remains unchanged in the 

Ashkin- Teller model. 
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INTRODUCTION 



Two dimensional statistical mechanics systems play a central role in our present understanding of phase transitions 
and critical phenomena. Outstanding members of this family of theories are the Ising model, the Ashkin- Teller 1 
and the eight-vertex- 2 - models. These models are useful to shed light on a variety of phenomena, in both classical 
and quantum physics, ranging from biological applications 3 - to the theory of cuprate superconductors^. Moreover, 
important advances in material science, accomplished over the last decades, have developed the ability to grow 
and experimentally explore ultrathin ferromagnetic films 5 , giving the opportunity to test some of the theoretical 
predictions. One of the fundamental questions concerning these essentially 2D materials is the role of defects and 
impurities in the critical properties of magnetic systems. Apart from academic interest, a detailed knowledge on the 
influence of defects on physical properties is always useful on general grounds, since all real materials are, to some 
extent, defected. In some cases of applied interest, such as ultrahigh-density magnetic recording media, it has been 
shown that linear defects can be used to efficiently control domain wall pinning, thus stabilizing the large area domain 
structure of ultrathin films^. Linear charge defects may also appear in graphene grown by chemical vapor deposition 
on Ni surface^. 

On the theoretical side, very little is known exactly about the behavior of planar systems in the presence of line 
defects^. For the simple square Ising lattice with an altered row (Bariev's model^) it has been shown that the scaling 
index of the magnetization varies continuously with the defect strength^^i -, whereas the critical exponent of the energy 
t-H ■ density at the defect line remains unchanged 11-13 . Taking this model as working bench, much insight was obtained 
I/"") ' about the origin of nonuniversal critical behavior™. More recently, by using path-integrals within the continuous 
^v^j . formulation of Ashkin- Teller and Baxter models, it was shown that the magnetic exponent depends on the strength 
7-H ' of the defect in exactly the same way as in Bariev's models. 

From another perspective, due to the well-known connection between the classical 2D Ising model and a quantum 
field theory of Dirac fermions in 1 + 1 dimensions, the study of defects as perturbations of conformal field theories 
has led to very important results in the area of integrable quantum field theories^—. This line of research was later 
focused on the problem of conductance in quantum wires^. The analysis of more mathematical aspects concerning 
the role of impurities and defects in the renormalization group flows of conformal models is currently under intense 
investigatio n 20 ' 21 . Very recently, the entanglement between two pieces of a quantum chain was analyzed by exploiting 
the connection with an Ising model with a defect line^ 2 -. 

All these advances were achieved for the case of homogeneous defects, i.e. when the defect strength is constant 
along the altered line. The case of non uniform couplings has been analyzed in the context of extended defects at 
surfaces 2 ^ - — and in the bul k 26 i 27 , displaying a rich variety of behaviors in the local critical properties. 

In this work we consider a narrow inhomogeneous defect and study the spin-spin correlator on the altered line. 
In other words, we analyze the extension of Bariev's model to the case in which the strength of the line defect is a 
function of the position on the column with modified couplings. Then, our result for the critical spin-spin correlator 
is a generalization of the result first obtained in Ref.[l3] for a uniform line defect. By using a path-integral approach 
in the continuum limit, we have obtained a formula that gives the spin-spin correlation function as a functional 
of an arbitrary defect distribution. This allows to explore the effect of different types of specific alterations in a 
straightforward way. We have also shown that the results remain valid for the Ashkin- Teller model, i.e. we found 
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that in these altered systems the non-scaling behavior of magnetic correlations on the inhomogeneous defect coincides 
with the one obtained in the Ising case. 

The paper is organized as follows. In Section II we explain our computational procedure for the well-known defect- 
free Ising model. In Section III we show how to extend the method when a line of altered couplings is included in the 
system. We emphasize how the case of inhomogeneous defect strength can be naturally considered with our technique. 
In Section IV we illustrate the use of our result showing the predictions for two specific defect functions. In Section V 
we extend the procedure to the more complex Ashkin- Teller and Baxter models. Finally in Section VI we summarize 
our findings and present our conclusions. 



II. THE METHOD: DEFECT-FREE CASE 



For completeness and illustrative purposes, we start by describing the computational procedure for the homogeneous 
defect-free case. The Hamiltonian of the original lattice model is given by 



H = -Y, J Wi (!) 



<v> 

where < ij > means that the sum runs over nearest neighbors of a square lattice (a = ±1). 

As shown in Ref. |28| the scaling regime of the 2D IM can be described in the continuum limit in terms of a model 
of Majorana fermions with Lagrangian density: 

C[a] = ai$a (2) 

where a represents a Majorana spinor with components ai,a- Let us recall that this components are connected to 

fermion annihilation and creation operators c r and c\. attached to site r (c r = e (cti (r) + ia^fV)). It is also useful 

for later convenience to define the energy-density as e Q = a\oi2- The symbol $ stands for 7„c? v , with j v the usual 
Euclidean Dirac matrices (y = 0, 1 associated to space directions). 

Similar manipulations, based on the Jordan- Wigner transformation™, allow to write the on-line spin-spin correlation 
function in the form— ^ 

< a(0)a(R) >isin 9 =< exp (tt / dxe a (x)) > (3) 

Jo 

where the vacuum expectation value is an anticommuting path-integral to be evaluated with the continuum action 
S = J d 2 x£, with an integration measure T>a. The explicit computation of ((3]) can be performed either in terms of 
the Majorana a-fields or in terms of Dirac fermions 30 built through the doubling technique 31 , yielding the well-known 
result for the Ising correlator. We start by squaring 

f R 

< ct(0)ct(R) >j smg = < exp (tt / dx (e a (x) + e a >(x))) > (4) 

Jo 

where the vacuum expectation value must now be computed with respect to an Euclidean action with Lagrangian 
density C[a, a'] — C[a] + C[a'], a' being the replicated fermion fields. Following Ref.[3fjJ we can build Dirac fermions 
^ as 

^ = a + ia'. (5) 
In terms of these new fields we can write the Lagrangian density C[a, a'] in the form 

= (6) 
where 75 = 17071 . On the other hand equation (j4]) can be expressed as 

< a(Q)a(R) >j sing =< exp (tt J d 2 x^£^ >, (7) 

where now the path integral integration measure in the right hand side is expressed in terms of the fields VP, and A v 
is an auxiliary vector field with components: 

A (xo,x 1 ) = -5(x )6[x 1 )6(R-x 1 ), A 1 (x Q ,x 1 ) = 0. (8) 
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Gathering the above results we can write: 



where 



< a(0)a(R) >j smg = f|^J, (9) 



Z[g] = J WW cxp(- J d 2 x(C[^} + g^A^)). (10) 

The continuum limit of the squared two-point spin correlation function is exactly expressed in terms of the vacuum 
to vacuum functional of a quantum field theory describing a Dirac fermion interacting with a classical background 
A v . Now we make the following change of path-integral variables in the numerator of equation with chiral and 
gauge parameters $ and rj, respectively: 

V]> = e - ?r (75*-i»?) £ ^ _ £ g-TT (75 

The integration measures T>*S> and 2?£ are related through the so called Fujikawa Jacobian J, V^'D'^' = 
</[$, 77] 2?£2?£. If the parameters of the transformation are related to the previously introduced vector field A v in 
the form 

A v = e vp d p <S> + d v r) (12) 
one easily gets Z[g = it] = J Z[g = 0], which leads to 

< a(0)a(R) > 2 Isin = J(R), (13) 



As explained in Ref. 32j, the Jacobian J(R) must be computed with a gauge- invariant regularization prescription 
in order to avoid an unphysical linear divergence. Following this procedure one finds that J depends on the $-field 
only as 

J(R) =exp-vr/2 J d 2 x d v <f>(x, R) d v $(x, R). (14) 

The explicit form of 3>(x, R) is determined by combining Eqs. © and (fl"2"|) which gives the following partial 
differential equation for $: 

n<S>{x ,x l ,R) = -5[x )^[0{x 1 )6{R-x 1 )] (15) 
ax 1 

where □ = Oq+O 2 . The solution of this equation is easily obtained by using the Green's function of the D'Alembertian: 
Gq(zo, zi) — j- ln(zg + z\ + a 2 ), with a an ultraviolet cutoff related to the original lattice spacing. Replacing in (fT4|) 
and considering the limit R » a we find the well-known result < a(0)a(R) >i s ing— (a/R) 1 / 4 . 



III. INHOMOGENEOUS LINE DEFECT 



Now we include a line defect in the original Ising lattice. To be specific we consider the so called chain defect (here 
we employ the terminology of Ref. [gj], which corresponds to Bariev's second type defect, in which bonds along the 
same column are replaced: Ji — > J' 2 )- In previous studies the altered coupling J' 2 was taken as a constant. From now 
on we allow J' 2 to vary from site to site, i.e. we make J 2 — > J 2 {x\). 

We will study the two-spin correlation function in the column of altered bonds (xq = 0)^. It is known that the 
continuous version of the classical model is modified, due to the defect, by the addition in equation ([2]) of a term 
2tt)i(xi) S(xq) e a (x), with /x = J 2 {x\) — J2 (see for instance 1131 ) . By carefully examining the fermionic representation 
of er-spin operators on the lattice, following the lines of Ref. [30(, one also finds that in the continuum limit each spin 
operator on the defect line picks up a similar /i-dependent factor, in such a way that the correlator for the defective 
model is given by a simple modification of equation (|3|): 

< a(0)a(R) >i nhom =< exp (tt I dxi (1 + 4/i(xi)) e a (x{)) >n ■ (16) 
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It is evident that the squared correlator can be written again as in equation ([7]). The presence of the inhomogeneous 
defect manifests in the form of the A^-field which is now rescaled by a factor (1 + 4/i(xi)). Thus equation (J5J becomes 

A (x ,xi) = -(l + 4^(xi))5(x o )0(x 1 )0( J R-xi) (17) 



Ai(a;o,xi)=0. (18) 

The implementation of the change of variables given by (| 1 1 1) and (TT21 leads to the generalization of (fT3")) : 

< a{0)a{R) > 2 mhom = Jinhom(R)- (19) 

Formally Jinhom(R) is still given by (|14p . but the effects coming from the nonuniformity of the defect strength enters 
the game through the function <fr(xo, x%, R), which now obeys a non trivial differential equation depending on /i(xi): 

n*(x ,xuR) = -S(x ) + 4/i(x x )) e{x x )6{R - Xl )\. (20) 

ax i 

The formal solution of this equation is 



®(x ,Xi,R) = i l n ^+f+^- fl ) 2 + 
+ 1 j R dx[ »(x[) £r In [xl + ( Xl - + a 2 ]. (21) 

Replacing in the corresponding expression for Ji n hom{R) we obtain 



< a(0)a(R) > mhom = {^ w )k+^^ 1 e p ( R \ (22) 

where 



F(R) = - f R dx u(x) — Tin (° 2 +(*-«) a ) (1+< "- W) 1 - 

- 1 So So R dxdy{\ + 4n(x)) i»(y)£[ln (a 2 + (x - y) 2 )] . (23) 

In the above integrals we have dropped the subindex 1 in the integration variables, in order to simplify the notation 
(xi — > x and y\ — > y). It is easy to check that in the special case fi(x) — > /z = constant, one obtains 

<a(0)o(R)>^(^f A °, (24) 

with A CT = |(1 + 4/i) 2 , which is the well-known result first obtained by McCoy and PerkiS. 

Formulae (|22[) and (|23l) constitute the main formal result of this paper. They give the critical spin-spin correlation 
on the altered line of an Ising model, as a functional of an arbitrarily varying defect strength. In the next section we 
will show some specific predictions for definite defect distributions. 



IV. APPLICATION TO SOME SPECIFIC DEFECTS 



Let us now consider some specific defect-functions for which F(R) can be analytically evaluated. We start with the 
following defect distribution 

V( x ) = p°7TT1 — FITS' ( 25 ) 
(1+ I x | /b) 

where b is a characteristic length scale. This function is similar to the one considered by Bariev in his study of 
horizontal large scale inhomogeneities^. Passing to dimensionless variables r = R/a and (3 = b/a, and considering 
weak defect strengths (/xo <C 1), for R S> a and 6>owe obtain 
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1 TO(2ff + r) / A. r-(2/3 + r) , f ^ 

We see that the magnetic correlation exhibits non-scalin g be havior, as expected for a local inhomogeneity. This 
result is in qualitative agreement with the analysis of Ref. [261 ] . However we should stress that we are considering a 
different situation here. Indeed, the present case corresponds to a standard 2D Ising model in which just one column 
(xq = 0) is altered in a non uniform way, whereas in Ref. [26j the couplings along columns are kept constant, while 
the couplings along all rows are modified in a non uniform fashion. In Figure [T] we compare the decays of correlations 
for constant defect (solid line), non constant defect with decay law ([25]) (pointed line) and the universal defect- free 
behavior (dashed line). In agreement with physical intuition the correlation decays monotonically with distance, in 
an intermediate way, faster than the defect- free case and slower than the case in which the defect strength is constant. 
Another expected feature, well reproduced by our solution, concerns the behavior with j3 — b/a: for increasing /3 the 
non-scaling decay becomes faster, being undistinguishable from the uniform case for large enough j3. 



<a(0)cr(r) > 

2r 
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FIG. 1: Spin-spin correlation as a function of r for constant (solid line) and variable (pointed line) defect strength. We set 
P = 10. The dashed line indicates the defect-free universal behavior r -1 / 4 . 



Let us now study a different function n(x), which represents a non-monotonic alteration of the line Xq = 0. For 
simplicity we consider a sequence of N slabs of heights fi oi (i — 1, ...N). Each slab starts at x = aci and ends at 
x = adi\ 

N 

fj,(x) — /ioi 0(x — aci) 0{adi — x), (27) 
i=i 

where 6(x) is Heaviside's function. Evaluating F(R) and replacing in (j22"j) . in the weak coupling regime (/zoi <C 1) 
and for r, c.i , ^ > 1 we obtain 

< <KowD ><—= C-rfl iC-fij- i^)]-*-™—/' [( ^(^£±i )r *-«' 2 (28) 

i— 1 

In Figure [2] we display the result given by the above formula for the simplest case: one slab or "barrier" starting at 
x/a = c = 10 and ending at x/a = d = 50. For r < c the critical two-spin correlation coincides with the standard, non 
defected correlation. In the presence of the defect, for c < r < d, it exhibits a faster decay. The correlation reaches 
a local minimum at r = d, and then it starts growing, approaching again the universal behavior corresponding to the 
magnetic critical index 1/8, asymptotically. In Figure [31 taking into account that ([2"8")l is valid for both positive and 
negative values of /zq, we show the critical correlation for a defect which is oscillatory along certain portion of the 
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line xo = 0, a sequence of five alternated slabs (/io = 0.1) and wells (fio = —0.1). As before, the spin-spin function 
coincides with the non defected one, for small distances (r < c\). For c\ < r < c? 5 there is an oscillatory behavior 
around the universal curve r -1 / 4 . For large distances the correlation tends to the universal decay. 



<a(0) a(r) > 




FIG. 2: Spin-spin correlation as a function of r for a line defect given by a slab starting at x/a = 10 and ending at x/a = 50, 
for /io = 0.1 (pointed line). The dashed line indicates the defect-free universal behavior r _1//4 . The solid line corresponds to a 
uniform defect with /io =0.1. 



<a(0) a(r) > 
1 
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FIG. 3: Spin-spin correlation as a function of r for a line defect given by an oscillatory function (a sequence of 5 slabs and 
wells with equal heights (depths) and widths (| no |= 0.1), starting at x/a — 10 and ending at x/a — 119) (pointed line). The 
dashed line indicates the defect-free universal behavior r -1 ^ 4 . The solid line corresponds to a uniform defect with /io = 0.1. 



V. EXTENSION TO THE ASHKIN-TELLER MODEL 



In this Section we show how to extend the treatment of inhomogeneous linear impurities to the Ashkin- Teller system 
consisting of two Ising subsystems with spin variables <jj and tj coupled by a quartic interaction 1 ™. The corresponding 
lattice Hamiltonian reads 

H = - ^ [ J2 (ViC] + TiTj) + J 4 O-jCTj-TiTj] (29) 

<ij> 

where < ij > means that the sum runs over nearest neighbors of a square lattice (a, t, = ±1). As it is well known, in 
the vicinity of the critical point this model can be described in the continuum limit in terms of two Majorana fermions 
interacting via their energy-densities: 



L[a, j3\ = ai$a + (3i0/3 — A e Q ep 



(30) 
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where a and f} are the Majorana spinors with components otip, f}\ t 2 respectively. e Q = u\U2 and ep = ftifa are the 
corresponding energy-densities. The coupling constant A is proportional to J4/J2. Let us now include, as before, a 
linear defect affecting one of the original Ising lattices, say the one with spins a. If this impurity is placed at column 
xq = 0, in the continuum limit we have to add to £ a term 2ivfi(xi) S(xq) £ a {x), with fi = J' 2 — J%. As shown in 
Ref. Il5l in order to compute the spin-spin correlator on the altered line is still possible to use the doubling technique 
depicted in Section II. However, in spite of the formal analogy, the situation is much more complex here. First of all, 
since we have two sets of spins, we have to introduce two Dirac fields: 4" = a + io! and \ = f} + if}' . We then obtain 



< a(0)a(R) > 2 AT =< exp (n J d 2 x *y4*^ > A 



(31) 



Here the background field A v is given by (|17[) and the vacuum expectation value must be computed with respect to 
an Euclidean action with Lagrangian density £[\&,x]: 

- ^[X75X*75* + /m( X T 7ix)^™(* T 7i*)], (32) 

where 75 = 17071 and \I/ T , \ T are the transposed spinors. 

The implementation of the change of variables given by pT|) and (TT2"j) leads to 

< a(0)a(R) > 2 AT =< a(Q)a(R) >j nhom F(A, R, n) (33) 

where < a(0)a(R) >i n hom is the defected Ising correlator given in (|22|) and 

F(X, R, n) = Af(\) < exp[^(C, x) + S v ((, X )] >o (34) 

where <>o means vacuum expectation value with respect to the model of free x an d C fermions. A/"(A) is a normal- 
ization constant independent of R. Since the analysis of the dependence of F(X,fi,R) on R is more easily done in 
momentum space, we have Fourier transformed S$(f, x, ji) and S n (^, X ,fJ>) in the above equation: 



s*(£x,/*) = ^ I f[-^[x(pih 5 x(P2)C(P3h 5 G(p 7 R,fi)C(p 4 ), 



(35) 



with G(P, R, fi) being a diagonal 2x2 matrix given by 



c<^>=( 9+( 7' M \^ *,,))' < 36 > 

where g± (P, R, fx) = ± J d 2 xe lPx e T2^<^(x,^,R) anc j p = p^^ + p 2 -\- p a -\- p 4 . A similar expression is obtained for S v with 
G{P, R) replaced by 

( m T%p,°,«,,o). ^ 

with h{P,R,n) = J d 2 xe iPx e 2 ^ x ' R ^. The explicit functional forms of <$>(x,R,(i) and r](x,R,fi) can be determined 
following the same steps depicted in previous Sections, yielding 



(x ,x 1 ,R^) = ^ln( ^^f ) + 



+¥Jo ax 1 n{x 1 ) (x i +{xi _ x , i)2+a2) 



(38) 



and 



v(^ 1 ,R, ( ,) = ^^dy T ^^ ) 



(39) 
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Then, g(P,R,fi) becomes 



, o td ixn I \/ Xn +a 2 arctan ( —n — ^ - — ) i4xo f^dy — n — rJ-rr^ tt 

h{P,R,n)= d 2 xe lPx e Vo V^ + ^+<^-W e 0J ° "(.g+^+^-i) 11 ) . (41) 



and h(P, R) 



Since any possible dependence on R of the function F(X, R, fi) comes from g±(P, R, /x) and h{P, R, fi), our problem 
is reduced to the analysis of these integrals. Let us first introduce a cutoff L, which can be interpreted as the size 
of the system, in order to avoid infrared divergencies (the thermodynamic limit will be recovered at the end of the 
computation by setting L —} 00). In terms of the dimensionless variable u p = x p /L, (p — 0, 1) we obtain 



g±(P,R,fj.) = lim ±L 2 / d 2 u 

L ^°° J\u \<l 



\u p \<l 

x r *PLu ( < + a 2 /L 2 + (ui - (R/L)) 2 ±1/2 Ti /^vKv) ^^bl)? 
1 u 2 + ayL2+u 2 > 

= ±(2ir) 2 5 2 (P) (42) 

and a similar result for h(P, R). Then, in the thermodynamic limit (a <C R <C L) F(X, R, fi) becomes independent of 
R and the critical behavior coincides with the one of the 2D Ising model in presence of an arbitrary inhomogeneous 
defect. 



VI. SUMMARY AND CONCLUSIONS 



We have considered the critical behavior of the two-spin correlation function in the continuum, field-theory version 
of the 2D Ising model with a line defect placed at the column xo = 0. In contrast to previous studies, here we have 
taken into account possible variations of the defect strength with the position on the line. Our main result (Eqs. ([22j) 
and (|23|l ) provides an analytical expression for the critical spin-spin correlation as a functional of an arbitrary defect 
distribution. From this one can explore the effect of different types of non uniform impurity distributions on the 
magnetization. In particular our finding can be used to analyze, within the critical regime, the transit from scaling 
to non scaling behavior. As examples, in order to illustrate the approach and check its validity, we have discussed 
two special cases: a defect strength decaying monotonously with distance from a given point, and a sequence of slabs. 
Finally, we extended the analysis to a nonhomogeneous line defect placed at one column of an Ashkin- Teller system, 
showing that the spin correlator on the altered line decays, in the thermodynamic limit, in the same way as in the 
Ising model. 
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